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ABSTRACT 

We report high-energy spectral data of the prompt emission of GRB970228 and 
its X-ray afterglow. We establish that the nature of the X-ray afterglow emission is 
non-thermal and similar to the later portion of GRB970228. Our data can be used 
to discriminate different emission models of GRB afterglows. While cooling of excited 
compact objects can be ruled out, fireball models are constrained in the physics of the 
radiation emission processes and their evolution. 



Subject headings: gamma rays: bursts — gamma rays: observations — stars: neutron 
— X-rays: general 
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1. Introduction 



Recent observations of the GRB970228 afterglow flCosta et al. 1997a ) by the BeppoSAX 



satellite (Boella et al. 1997) have shown that it may be actually possible to identify the ultimate 
source and nature of emission of gamma-ray bursts (GRBs). Eight hours after the burst, we 
detected a new transient X-ray source, 1SAX J0501. 7+1146, by re-pointing the BeppoSAX 
Narrow Field Instruments (NFIs) at the GRB location. Various arguments render compelling 
the interpretation of this source in terms of a fading X-ray afterglow from GRB970228: (a) the 
backward extrapolation of the source light curve to the time of the burst is consistent with the 
average X-ray flux detected during the late phase of the burst ( Costa et al. 1997b|) ; (b) its celestial 
position is consistent with the burst position determined by the Interplanetary Network ( |Hurley| 
et al. 1997| ); (c) a fading optical transient consistent with both the Wide Field Camera (WFC) 
and NFIs positions was also discovered ( [van Paradijs et al. 1997 , Guarnieri et al. 1997j ); (d) an 
observation of the GRB error box with ROSAT provided the position of the X-ray transient 
source within a 10" radius. Its centroid coincides, within 2", with the optical transient position 
( Frontera et al. 1997b , Frontera et al. 1997cj ). The optical transient appears to be embedded in an 
apparently constant extended source ftSahu et al. 1997|) . If it is a galaxy, those models that place 
GRBs at cosmological distances would be favoured. 



In another paper flCosta et al. 1997b ) we reported the decay of the light curve of the source, 



which is described by a power law function (oc t a with index a = 1.33^q that is consistent 



with that expected in simple relativistically expanding fireball models ( |Wijers et al. 1997 ). In this 
paper we report on the spectral evolution of GRB970228 and its afterglow. 



2. Observations 



GRB970228 was detected by the Gamma-Ray Burst Monitor (GRBM, 40-700 keV) ( |Frontera 



et al. 1997a| ) and WFC No. 1 (1.5-26.1 keV) flJager et al.T997j ) aboard BeppoSAX on February 



28.123620 UT QCosta et al. 1997a| ). Its position was determined to be within an error radius of 
3' (3cr) centered at a 2 ooo = 05 h 01 m 57 s , <5 2 ooo = H°46 / 24". Eight hours after the initial detection, 
the Narrow Field Instruments on-board BeppoSAX were pointed at the burst location for a first 
target of opportunity (TOOl) observation from February 28.4681 to February 28.8330 UT. A new 
X-ray source, 1SAX J0501. 7+1146, was detected (|Costa et al. 1997b|) in the GRB error box by 



both the Low Energy (0.1-10 keV, 8725 s exposure time) and Medium Energy (2-10 keV, 13444 s) 
Concentrators Spectrometers (LECS and MECS). The same field was again observed about three 
days later for a second observation (T002) from March 3.7345 to March 4.1174 (16270 s MECS 
and 9510 s LECS). During this observation, the 2-10 keV source flux had decreased by about 
a factor 20. The main arguments for the association of this source with the X-ray afterglow of 
GRB970228 have been discussed in Section 1. 



The data available from GRBM for spectral analysis include two 1 s ratemeters (40-700 keV 
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and >100 keV) and 128 s count spectra (40-700 keV, 225 channels). The energy resolution of the 
GRBM unit 1, co-aligned with the WFC No. 1, is 20% at 280 keV flAmati et al. 1997Q . WFCs 
(energy resolution ~ 20% at 6 keV) were operated in normal mode with 31 channels in 1.5-26 keV 
and 0.5 ms time resolution QJager et al. 1997 ). The burst direction was offset by 14° with respect 



to the WFC axis. The effective area of the GRBM unit 1, co-aligned with the WFC No. 1, is 
~ 420 cm 2 in the 40-700 keV band and is ~ 500 cm 2 at 300 keV. The corresponding effective 
area of WFC No. 1 averaged in the 2-26 keV energy band is 118 cm 2 . The background level 
in the WFC and GRBM energy bands was variable during the burst ( [Feroci et al. 1997| ). This 
variation was apparent up to 100 keV, while above 100 keV the background was stable. For the 
WFC spectra and light curves, the background level was estimated using an equivalent section 
of the detector area not illuminated by the burst. In the case of GRBM, background during the 
burst was estimated by interpolation, using a quadratic function that fit the 150 s count rate 
data before the burst and the 100 s data after its end. The presence of a variable background 
prevented the derivation of useful upper limits to the afterglow emission after the burst. During 
TOOl, the X-ray source associated with the burst afterglow was offset by about 10' with respect 
to the telescope axes. Since the image centroid, for MECS 1 and 2, coincided with of a strong-back 
that supports the detector window and absorbs most photons up to about 6 keV, we used for 
the spectral analysis only MECS 3 and LECS. The LECS and MECS spectra for TOOl were 
extracted from a 4' radius region around the centroid of the image source. For T002, the source 
was on-axis and the spectra from LECS and MECS were extracted from a 2' radius region. The 
MECS spectra were equalized and added together. The background spectrum was estimated from 
long observations (200 ks for MECS and 100 ks for LECS) of blank sky fields, using an equivalent 
detector region of the source image with a similar offset. 

Figure [l] shows the time profiles of GRB970228 in different energy bands after the background 
subtraction. In the 7-ray energy band (40-700 keV), the burst is characterized by an initial 5 s 
strong pulse followed, after 30 s, by a set of three additional pulses of decreasing intensity. In the 
X-ray band (1.5-10 keV) the first pulse starts within one second of the 7-ray pulse, achieves the 
peak flux later and has a duration three times longer. The entire burst duration is ~ 80 s. The 
7-ray (40-700 keV) fluence of the entire burst is (1.1 ± 0.1) x 10~ 5 erg cm~ 2 , which is a factor 
5 ± 1 higher than the corresponding value in the 2-10 keV band. (This ratio is based on the 
more accurate analysis of the WFC data performed for this paper and corrects the approximate 
value given previously by Costa et al. 1997c) . For comparison, the 40-700 keV fluence integrated 
over the first pulse is (6.1 ± 0.2) x 10~ 6 erg cm~ 2 , which is a factor of 8.5 ± 1.3 higher than the 
corresponding fluence in 2-10 keV. The 7-ray peak flux is (3.7 ± 0.1) x 10 -6 erg cm -2 s , while 
the corresponding 2-10 keV flux is (1.4 ±0.1) x 10~ 7 erg cm" 2 s , with a ratio between 7-ray and 
X-ray peak fluxes of 26 ± 2. In the BATSE energy band (50-300 keV) the fluence of GRB970228 
is 6.1xl0 -6 erg cm -2 . 
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3. Spectral analysis and results 

We divided the GRB time profile into seven temporal sections, and performed a separate 
spectral analysis on the average spectrum of each section. The duration of the sections were chosen 
on the basis of the statistical quality of the data. In particular, we obtained two spectra (5 and 9 s 
duration, respectively) for the 1 st pulse (A and B in Fig. [l]), four 3 s spectra for the 2 nd pulse (C, 
D, E and F), and only one spectrum for the 3 rd and 4 th pulses together (G in Fig. p]). All WFC 
spectra are well fit with power laws {Kx ■ E ax ), while black-body spectra do not fit the data. We 
assumed power laws (Kr ■ E ar ) to deconvolve spectral data in the GRBM band. In Fig. || we show 
the logarithmic power per photon energy decade (the vF v spectrum) for six of the time intervals 
defined above. The spectral evolution of the 1 st GRB pulse and its diversity from the subsequent 
pulses is evident. The initial part of the l si pulse (Fig. ^A) is quite hard, with noticeable curvature 
of the vF v spectrum. The peak energy, that is not detectable in the WFC/GRBM energy range, 
is likely above 700 keV. The decay part of the I s * pulse shows substantial evolution, with a clearly 
detectable broad curvature. By fitting the photon spectrum with a smoothly broken power law 
QBand et al. 1993| ), we found a peak energy at E p = 35 ± 18 keV. This peak energy behaviour is 



typical of GRB spectral evolution ( Ford et al. 1995 ). In Fig. |3| we show the time evolution from 



the burst onset of the GRBM and WFC power law spectral indices «r and ax- For the 1 st pulse, 
the temporal behaviour of a?r is well described by a linear law (ar(t) = a + bt, \ 2 = 1-28, 3 dof) 
with a = —0.9 ± 0.2 and b = —0.31 ± 0.05, and t in seconds. The power law photon indices ctr and 
ax do not evolve in the same way as a consequence of the curvature of the I s * pulse spectra. We 
also note the interesting spectral softening during the decaying part of the pulse. 

The spectral behaviour of the 2 nd pulse is significantly different from that of the 1 st one. The 
spectral data can be fit with a single power law over the 1.5-700 keV range. During the first 3 s of 
the 2 nd pulse the spectrum is significantly harder than during the last part of the l si pulse. The 
photon index or of the 2 nd pulse vs. time from pulse onset shows some evidence of evolution (it 
can be fit with a linear law as above with a = —1.6 ± 0.2 and b = —0.12 ± 0.04), but it is also 
consistent with a constant value. Instead ax is definitely constant during the 2 nd pulse. The 
overall WFC/GRBM photon index of the 2 nd pulse is initially -1.8 ± 0.1 and then -1.9 ± 0.1. 
Interestingly, there is no strong evidence for spectral curvature as shown in Fig. ^, probably 
indicating that the peak energy E p passed through the WFC/GRBM energy band during the first 
part of the burst and, as a consequence, later it is below ~ 2 keV. Low counting statistics prevents 
us from deriving good quality spectra for the last section of the 2 nd pulse (section F) and for the 
3 rd and 4 th pulses, separately. For section F (not shown in Fig. ^), the X-ray data were fit by a 
power-law model with ax = —1-5 ± 0.4, while the 7-ray data provided only a 2a upper limit of 
the ar spectral index of —0.6. The average photon spectrum of the sum of the 3 rd and 4 th is well 
fit with a power-law model of ar = — 1.4^q 3 and ax = —1.6 ± 0.1. It appears that the spectrum 
is again harder than that at the end of the 2 nd pulse and more similar to the one at the beginning 
of the 2 nd pulse. In any case, it is consistent with a single power-law. The vF v spectrum of the 
sum of the 3 rd plus 4 th pulses is shown in Fig. |2|G. 
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Remarkably, we find that the power-law spectral model continues to hold also during the 
X-ray afterglow. The photon spectrum of 1SAX J0501. 7+1146 in the 0.1-10 keV energy band 
as measured during the TOOl observation is shown in Fig. |I[ It was obtained from the source 
count spectra of LECS and MECS after deconvolution for the instrument responses. Spectral 
fits were performed with the XSPEC 9.0 package. We fit various functions to the source count 
rate spectrum (see Table ||). The best fit model, with a x 2 = 7.2 (10 dof) is a power-law 
of photon index a = 2.06 ± 0.24 and photoelectric absorption N# = 3.5^2 3 x 10 21 cm~ 2 . A 
power-law model with N# = does not fit the data nor does an absorbed black-body model. 
The probability that the difference in the x 2 between these models and the absorbed power-law 
is due to chance is 1.2 x 10~ 3 and 7.2 x 10 -4 , for the power-law with N# = and for the 
absorbed black-body, respectively. Our determination of N# is consistent with the total Galactic 
photoelectric absorption (~ 1.6 x 10 21 cm -2 ) expected along the line of sight to the GRB error 
box. We can exclude a nearby source with zero equivalent N# (e.g., near the Solar system) at 
about 3a statistical significance. In Fig. [|, we also show the best fit model spectrum of the X-ray 
afterglow. During the T002 observation, the X-ray source was too weak to derive a high-quality 
spectrum. The MECS detected a 2-10 keV flux of (1.5 ±0.5) x 10~ 12 erg cm -2 s _1 from the source, 
while the LECS did not detect any significant flux (3a upper limit of 0.4 x 10 -12 erg cm -2 s _1 ). 
The hardness ratio HRi between the count rate in the hard X-ray band (2-10 keV) and that in 
the soft band (0.1-2 keV) is 2.5 ± 0.8 for TOOl and >2.5 for T002 (3a lower limit). However 
the hardness ratio HR2 between two hard X-ray bands (3-5 keV and 1.5-3 keV) does not show 
statistically significant variations from the first to the second TOO observation. Thus we cannot 
infer a change in the spectral shape from TOOl to T002. 

4. Discussion and conclusions 

The non-thermal spectrum of the burst radiation and of the X-ray afterglow is the most 
important result of this analysis. A power law model was found to fit the spectral data better 
than a black body model. The maximum value of or (—0.9 ± 0.2) achieved at the onset of the 
I s * pulse as well as the corresponding ax during the first 5 s are consistent with the asymptotic 
spectral index below the peak energy E p (—0.67), as expected in synchrotron emission models 
(e.g., the Synchrotron Shock Model (SSM) ( [Tavani 1996[ )). The spectral index of the I s * pulse 
rapidly evolves to Qr = —2.3 ± 0.3 at the end of the 7-ray pulse. In the SSM framework the 
final value of ar during the 1 st pulse corresponds to an index of the non-thermal electron energy 
distribution function 5 = 3.7. During the first 3 s of the 2 nd pulse the spectrum is significantly 
harder than during the last part of the I s * pulse. This indicates significant re-energization or 
relaxation of the particle energy distribution function within the ~ 20 s between the 1 st and the 
2 nd pulse. It appears that the last three GRB pulses and the X-ray afterglow have a similar 
non-thermal spectrum and that this spectrum does not appear to change from the first to the 
second TOO. These results are not consistent with simple cooling models of excited compact 
objects. An analysis of Ginga data suggested the existence of black-body spectral components in 
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the 'precursor' or 'delayed' GRB emission in the 1-10 keV band (Yoshida et al. 1989, Murakami 
et al. 199l| ). In this burst we find no evidence for a black-body spectral component of temperature 
~ 1-2 keV or for a prominent soft X-ray component. Furthermore, there is no evidence of upturn 
in the soft X-ray intensity with respect to the higher energy spectrum (Liang et al. 1997). 



An evolving non-thermal spectrum for both the burst and afterglow emission is generally 
expected in relativistic expanding fireball models ( Meszaros and Rees 1997 , Tavani 1997). In these 
models, the physics and locations of the shocks associated with the fireballs heavily influence 
the photon emission mechanisms. Simple fireball models, in which only the forward blast wave 
radiates efficiently, predict an evolution of the peak energy as a function of the time t from the 



burst onset, E p oc t 3 / 2 ( Tavani 1997 , Wijers et al. 1997 ). By extrapolating our data for the 1 



I st 



pulse of GRB970228 we obtain an initial E p ~ 1 MeV. In this model, the overall evolution of the 
X-ray intensity is expected to evolve as oc t s with 5 = (3/2)(a + 1) and a an appropriate photon 
index ( Wijers et al. 1997 ). If we use our best value of «r for the decay part of the 2 nd pulse 
(—1.94 ± 0.13), we obtain 5 = (—1.4 ± 0.2), a value that is consistent with the observed decay of 
the GRB970228 afterglow flCosta et al. 1997b| ). However, if we use the value of ar determined 
at the end of the 1 st pulse (—2.3 ± 0.3), the resulting value of 5 (—1.95 ± 0.45) would not agree 
with our observations. The discontinuity in the 7-ray spectral index observed from the end of the 
I s * pulse to the beginning of the 2 nd pulse requires an interpretation. Given the continuity of the 
spectral index, starting from the 2 nd pulse, any physical relation between the X-ray component of 
the GRB and the afterglow emission most likely holds with the last set of hard GRB pulses, not 
with the I s * one. This suggests that the emission mechanism producing the X-ray afterglow might 
be already taking place after the 1 st pulse. 



This research is supported by the Italian Space Agency (ASI). KH is grateful to the US 1SAX 
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Table 1. Spectral fits of 1SAX J0501. 7+1146 in the energy range 0.1-10 keV 



Model A b a kT (keV) N H (xl0 21 cm- 2 ) X 2 dof 

Power Law plus abs c l-l^io 2.06 ± 0.24 ••• 3.5±|J 7.2 10 

PL with abs set to 0.52 ± 0.15 0.45 ± 0.20 • • • (frozen) 17.8 11 

Black Body plus abs c 0.032 ± 0.005 ••• 0.86 ±0.12 +1 18.7 10 



a Uncertainties at the 90% confidence level for a single parameter. 
b Photons cm -2 s _1 keV^ 1 at 1 keV. 

c Photoelectric absorption of a gas with cosmic abundance ( Morrison and McCammon 1983j ). 
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Fig. 1. — Light curve of GRB970228 in different energy bands, after background subtraction (see 
text). Dashed lines give the background reference. In the insets, the three minor pulses of the 
burst. The time sections on which the spectral analysis was performed are shown on the top. 
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Fig. 2. — vF v spectrum of the burst (u is the photon energy in keV and F v is the specific energy 
flux in keV cm~ 2 s _1 keV -1 ) for 6 of the 7 time sections in which we divided the burst time profile. 
The time intervals from the burst onset over which the spectrum has been accumulated are also 
shown. 
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Fig. 3. — Behavior of the X-ray and 7-ray power law photon indices with time from the burst 
onset. The photon index of the X-ray afterglow spectrum is also shown. 
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Fig. 4. — Deconvolved LECS and MECS spectra of the X-ray afterglow during TOOL The best 
fit power law model with low energy photoelectric absorption is superposed to the LECS (dashed 
line) and MECS (full line) data. 



